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ABSTRACT 

We describe the methodology and compute the illumination of geometrically 
thin accretion disks around black holes of arbitrary spin parameter a exposed to 
the radiation of a point-like, isotropic source at arbitrary height above the disk 
on its symmetry axis. We then provide analytic fitting formulae for the illumina- 
tion as a function of the source height h and the black hole angular momentum 
a. We find that for a source on the disk symmetry axis and h/M > 3, the main 
effect of the parameter a is allowing the disk to extend to smaller radii (approach- 
ing r/M — > 1 as a/M — > 1) and thus allow the illumination of regions of much 
higher rotational velocity and redshift. We also compute the illumination profiles 
for anisotropic emission associated with the motion of the source relative to the 
accretion disk and present the fractions of photons absorbed by the black hole, 
intercepted by the disk or escaping to infinity for both isotropic and anisotropic 
emission for a/M = and a/M = 0.99. As the anisotropy (of a source approach- 
ing the disk) increases the illumination profile reduces (approximately) to a single 
power-law, whose index, q, because of absorption of the beamed photons by the 
black hole, saturates to a value no higher than q > 3. Finally, we compute the 
fluorescence Fe line profiles associated with the specific illumination and compare 
them among various cases. 



Subject headings: accretion, accretion disks — black holes physics — neutron 
stars — stars: individual X-rays: reprocessing 
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Introduction 



The "standard" geometry of accreting active galactic nuclei (AGNs) (and for that matter 
also galactic black hole candidates) consists of an optically thick, geometrically thin accretion 
disk that extends to the radius of the innermost stable circular orbit (ISCO) surrounded by 
a hot (T ~ 10 8 — 10 9 K), X-ray emitting corona. Coronal X-rays are reprocessed on the 
surface of the much cooler disk (T ~ 10 4 — 10 6 K) to produce fluorescent Fe Ka emission line 
at 6.4 keV as well as a reflection component; the latter constitutes a broad peak near E ~ 30 
keV resulting from the reprocessing of the power law coronal emission on the colder disk by 
reducing its flux by photoelectric absorption at low energies and by down-Comptonization 
at high energies. 

The intrinsically narrow shape and the well known energy (6.4 keV) of the Ka line make 
it a unique diagnostic of the inner regions of accretion flows around black holes: An accretion 
disk such as those described above is expected to produce line profiles of unique shapes be- 
cause of the large Doppler mot ion and redshifts assoc iated with emission from the innermost 
radii of such disks (see however iTitarchuk et al.ll2003l . for an alternative view). Model fluores- 



cence Ka profiles ha ye been computed for disks around e ither Schwarzschild (|Fabian et al. 



198 91) or Kerr (e.g.. lLaorl Il99ll : Beckwith fe Pond 12004 ; iDovciak. Karas. fe Yaqoobl 12004 



Brenneman fe Reynoldsll2006l ) black holes, whi le more general treatment including the reflec- 
tion component has also been discussed (e.g.. iGeorge fe Fabianlll99ll : Matt. Perola. fe Piro 
199lh . 



The advent of observation has s ubsequently vindicated these considera tions with the 
discovery of broad iron Ka lines (e.g.. iTanaka et al.lll995l : iNandra et al.lll997l ). Of particular 
interest has been the Ka line associated with the AGN MCG -6-30-15, which, because of 
its very broad red wing, indicated emission by matter interior to the radius of the ISCO of 
a Schwarzschild black hole (r/M = 6 where M is the black hole ma ss), thereby implyin g 
emission by an accretion disk around a highly rotating Kerr black hole (jlwasawa et al.lll996l ). 
Observations of the same object by the much larger area of XMM-Newton incr eased the sig 



nificance of these observ ations and confirmed the presence of a very broad line (IWilms et al. 



200ll ; iFabian et al.ll2002l ). however the low sensitivity of the instrument above > 10 keV has 
made difficult the determination of the underlying contin uum and therefore the true width of 
the line. More recently, observations of MCG -6-30-15 by lMiniutti et al.l (120071 ) with Suzaku, 
with improved sensitivity at higher energies over XMM-Newton, were able to confirm that 
the Fe Ka line in this object is indeed as broad as it was originally claimed. 

Clearly, the profile of the Fe Ka emission depends, besides the geometric properties 
of the accretion disk, also on its kinematic properties and on its illumination by the X- 
ray continuum source. While the aforementioned treatments, as well as most of similar 
treatments in the literature assume axisymmetry and Keplerian motion for the plasma in the 
disk, these assumptions are not guaranteed by any fundamental principle. Non-axisymmetric 
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structures and/or non-negligible r adial velocity in th e disk can yield l i nes of different shapes 
as sh own by some authors (e.g., 



Ruszkowski 


2000 




Yu & Lu 


2000 





Fukumura & Tsuruta 



2004 ). Finally, even in axisymmetric geometries , different i llumination laws can yield line 



profiles of different shapes, as originally shown by lLaorl (1199 lh . For computational simplicity. 



Fabian et al. 


2002 


)• 


Laor ( 



illumination laws as well as with that of an infinite plane by a point source at height h in 
flat space, i.e., F(x) = h/[Aii(h 2 + x 2 ) 3//2 ], where x is the distance on the disk plane from the 
foot of the vertical line from the source to the plane. 



The sensitivity of the line profiles on this parameter (see lLaorl Il99ll . Fig. 2) has led 
us to take a closer look into this issue. This issue has been addressed to a certain extent 



by iReynolds fc Begelmanl (119971 ) who computed the illumination for a Schwarzschild black 
hole assuming an isotropic source at a given height, with emphasis on the contribution to 
the line emission by matter in free fall region (i.e., plunging region), interior to the ISCO at 
r/M = 6. The question we are set to study in this note is the illumination law of an accretion 
disk by isotropic and non-isotropic point sources at arbitrary values of the height h and for 
arbitrary values of the black hole spin parameter a. Our goal is to produce fitting formulae 
for the most generic cases which can be easily incorporated into existing numerical/software 
codes to produce more realistic line profiles for given source configurations. We place some 
emphasis on the issue of isotropy which changes in the presence of strong gravitational field. 
We also provide examples of anisotropic illumination by assuming that the source is moving 
instantaneously and that all anisotropy is due to the relativistic beaming of its photons. 

In §2 we outline the setup of the problem and the method followed for the case of a 
Schwarzschild black hole, while we defer discussion of the Kerr geometry to the Appendix. 
In §3 we present the results of our calculations along with their fitting formulae in graphic 
form for comparison with the direct calculation and provide a few examples of the resulting 
Fe Ka line profiles. Finally in §4 we conclude by summarizing and evaluating our results. 



2. Formulation of the Problem 

We consider an X-ray source located at a height h above the plane of an accretion disk 
and along its axis of symmetry as shown in Figure [H A black hole is located at the origin of 
the coordinate system, with the accretion disk extending to the ISCO. The latter depends 
on the spin parameter of the black hole, ranging from r/M = 6 for a Schwarzschild black 
hole (a/M — 0) to r/M = 1 as a/M — » 1. We also assume the X-ray emission to be isotropic 
in the local frame of the X-ray source. 



In this part of our treatment we follow the exposition given in 



Chandra sekhar 



(1983) 
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Fig. 1. — Schematic illustration of our model. An accretion disk is illuminated by a point-like 
X-ray source at height h on the symmetry axis. 



in formulating the problem of illumination and the method followed; we restrict ourselves 
to discussing the Schwarzschild case only, deferring the discussion of the Kerr geometry 
to the Appendix. We adopt the usual geometrized units in the Kerr geometry in which 
G = c = 1 where G and c are the gravitational constant and the speed of light, respectively. 
Thus, distance (length) is scaled by the black hole mass M throughout the paper. The 

M + VM 2 - a 2 



event horizon is then expressed as rn 
Schwarzschild case the radial equation of the photon orbit is 



2M as a/M -> 1. In the 



dr\ 2 9 L 2 ( 2M 



to be considered along with the angular momentum L and energy E integrals of motion 



r^ = L and 



2M 



d ± = E . 

dr 



(2) 



These, through the redefinition of the null affine parameter r to dr = Edr, read 



r 2 — = b and 
dr 



2M\ dt 
r ) df 



(3) 



where h = L/E is the impact parameter of the orbit, its unique degree of freedom. In the 
remaining of the paper we drop the tilde in the definition of r. 

On dividing equation §Q by r A (d6 / dr) 2 = L 2 one obtains the equation for the photon 
orbits 



1 dr 
r^d6 



1 

V 2 



R 2 (r) 



(4) 
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where R 2 (r) = r 2 /(l - 2M/r). 

As shown in Figure [T] we co nsider photons emitte d by the source at (local) polar angles 



ip defined by the expression (see IChandrasekharlll983l . p. 127) 



1 dr 1 dr 

COt ^ = rd9 = r(l-2M/ry/*d§ ' (5) 

where df = (1 — 2M/r)~ 1 / 2 dr is the element of proper length in the radial direction and 
r and 9 are the usual Boyer-Lindquist coordinates of the Schwarzschild geometry. We use 
the above equation (jSJ) for the definition of isotropic photon emission at the local frame of 
the source, and so we choose their distribution to be uniform in the angle ip (i.e., we set 
Aip = constant). Using the equation of the photon orbits [equation (j3J] and the definition 
of R(r), the above expression for the angle ip can be rewritten as 

R 2 (r) r 2 

cot 2 ip = —- 1 or b 2 = R 2 (r) sin 2 ^ = — sin 2 ?/' . (6) 

(r (1 — 2M/r) 

The above equation (j5J) provides the necessary relation between the photon impact parameter 
b and its local emission angle ip. 

The photon trajectories of a (locally) isotropic emission are determined by choosing 
the value of ip uniformly between and n and then using equation (jHJ) to determine the 
value of the photon's impact parameter b which, substituted in equation (TjJ, provides the 
associated photon orbits. However, because the photon orbit equation involves the square 
of the derivative dr/d8, we have found it more convenient to use the second-order equation 
for r supplemented with a value for the derivative of the radial coordinate, f (= dr/dr), 
to be used as an initial condition. Differentiating equation (PQ) with respect to r and using 
equation (jBj) we obtain the following expressions for the equation of motion and the initial 
condition on r 

3M\ „ b 2 b 2 ( 2M\ „ 
1 and r 2 = 1 - — — = 1 - — 1 = cos 2 ^ . (7) 



Cast in the above form, the photon equation of motion makes apparent the significance of 
the r/M = 3 surface for photon circular orbits in the Schwarzschild geometry; the second 
relation then indicates that, because of the single degree of freedom of the photon orbits, the 
value of r necessary for their specification is uniquely determined by the value of the impact 
parameter b or the angle if) through this relation. 

Our procedure for the determination of the orbits is the following: We consider a large 
number of orbits determined by the photon emission angle ip chosen uniformly between 
and 7T. Using equations ([6]) and (j7j) we obtain for each value of ip the values for b and r 
which are used to integrate the second-order photon equation of motion along with that for 
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the angular momentum, dd/dr = b/r 2 . The integration proceeds from an initial value 9 = 
(the location of the source on the symmetry axis) until the angle 9 reaches the value 9 = n/2 
(the position of the equatorial plane where the accretion disk lies) at which point we note 
the value of the radial coordinate r (now on the surface of the disk) at which the photon 
is located. Because we launch photons only on the poloidal plane from a source located on 
the symmetry axis of the source - black hole configuration, in order to take into account the 
photons emitted in the azimuthal direction, we weigh each photon intercepting the accretion 
disk by sim/>, the sine of the local polar photon emission angle. We collect all the (weighted) 
photons arriving at the disk plane (i.e., 9 = 7r/2) at a given radial interval, from the horizon 
rn to r/M = 2000 which is then converted to photons per unit area to produce the disk 
illumination. 



3. Results 

Following the procedure outlined above for a Schwarzschild black hole and for a Kerr 
black hole in the Appendix we have produced the illumination of accretion disks for a variety 
of values for the black hole spin parameter a and height of the source, h, expressed in units 
of the black hole mass M. Our results along with the fitting formulae appropriate for a 
most generic case, namely that of an isotropic point source on the symmetry axis, are given 
in the next subsection (§3.1). The effects of a locally anisotropic emission are discussed in 
the following subsection (§3.2), while we conclude this section by presenting the line profiles 
that result from the corresponding model illumination laws (§3.3). 



3.1. Illumination by an Isotropic X-ray Source 

We first present our results the case of an isotropic source on the axis of symme- 
try. In Figure [2] we show several of the calculated photon trajectories in the poloidal 
plane (r,9) for two cases: (a) h/M = 100 and (b) h/M = 3. In each case the source 
is isotropic locally and the photons are chosen in constant angular size intervals of size 
Ai/j = 8.7 x 10~ 2 radians = 5deg (corresponding to a total of iV = 36 orbits in the fig- 
ure). Because of the axial symmetry of the problem we only need show the trajectories of 
the photons emitted over < if) < tt. The value of black hole angular momentum used 
in both cases is a/M = 0.99. For an X-ray source at a height much larger than the black 
hole horizon (i.e., h/M ^> 1) as in case (a), the disk illumination is not very different from 
that of the Minkowski geometry given above [oc h/(h 2 + r 2 ) 3//2 ] (see e.g., Fig. [5^, also to be 
discussed later). On the other hand, as the source approaches the black hole, as in case (b) 
of Figure El the focusing influence of the geometry becomes significant; the light rays bend 
toward the black hole and enhance considerably the disk illumination in its vicinity; at the 
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same time, an increasing fraction of the photons are absorbed by the black hole with the 
fraction approaching unity as the source position approaches the horizon. It should be noted 
that in case (b) a significant fraction of the photons are also subjected to frame-dragging 
by the black hole angular momentum; this causes their rotation in the azimuthal direction 
before reaching the horizon (note the azimuthal trajectories partially seen just outside the 
horizon denoted by a dark hole). However, because of the axisymmetry of the configuration, 
the frame-dragging effects that influence only the azimuthal component of the trajectory are 
ignorable; the illumination depends on the photon trajectories on the poloidal plane. 

The main effects of the black hole angular momentum a, as far as our current consider- 
ations go, mainly lie in the reduction of the sizes of the horizon rn and, most importantly, 
that of ISCO, which allows the illumination of regions of the disk with sufficiently large red- 
shift, a necessary condition for producing the very broad Fe lines observed in MCG -6-30-15. 
While the photon orbits in the poloidal plane can change significantly with increasing a, the 
illumination law depends mainly on the height of the source above the disk h/M and it is 
largely independent of the value of a (this is true for h/M > 3; for smaller values of this 
parameter, the illumination laws do differ substantially). We have verified the above state- 
ment by direct orbit computation: Figure [3] exhibits the illumination profiles of accretion 
disks around a Schwarzschild (solid) and a Kerr black hole of a/M = 0.99 (dotted curve) 
respectively by a source at h/M = 3. Even for such a small height, the illumination patterns 
are almost identical, except for the fact that a high value of a allows smaller value of the 
horizon rn and of ISCO. Because of this the illumination profiles obtained for a Kerr black 
hole are also appropriate for a black hole with any value of a, provided that the illumination 
is restricted to the proper value of ISCO. 

By computing a sufficiently large number of orbits in a Kerr geometry with a/M = 0.99, 
similar to those appearing in Figure [2], and keeping track of their normalization (proportional 
to simp) we have calculated the radial distribution of the illumination by collecting the 
photons intercepting the disk plane within a given range in radius. The results for different 
values of the parameter h/M are shown in Figure HJ the values of this parameter used 
are: h/M = 100,50,30,10,6,5,4 and 3. The inner radius is set to r-^/M = 1.15 (for 
a/M = 0.99) and the outer radius is set to r out /M = 1000. In all cases, the illumination 
approaches asymptotically (i.e., for r 3> h) the functional form appropriate for Minkowski 
space [oc h/{r 2 + h 2 ) 3 ^ 2 ~ r~ 3 for r ^> h\\ this result is not unexpected as the photons that 
intercept the outer reaches of the disk have been emitted almost horizontally at the source, 
directions at which the influence of the focusing effects of the black hole are minimal for the 
values of h/M considered here. The normalization of the illumination in this regime of the 
radius r is very close to that of Minkowski space too except for the case h/M — 3 in which 
the gravitational effects "pull" the photon trajectories toward the black hole reducing the 
normalization below its flat space value (see also Fig. [5]i to be explained later). At values 
of r ~ h the profile "levels-off" to a slope which approaches that of flat space (oc r°). The 
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r sin0 rsinG 

Fig. 2. — Photon trajectories in the poloidal plane (r, #)-plane for (a) h/M = 100 and (b) 
h/M = 3 around a rotating black hole with a/M = 0.99. An equatorial thin disk is depicted 
by the thick gray line at 9 = n/2 along with a black hole at (0, 0). The gray curve around 
the black hole denotes the static limit (ergosphere) projected onto this poloidal plane. Scales 
are different for presentation purposes. 

effects of gravitational focusing begin to become noticeable for r/M < 7 (even for very large 
values of h/M), as shown in Figure Eh; their importance increases with decreasing values of 
h/M, as expected and shown in Figure EMSti, where the actual illumination is compared to 
that corresponding to the Minkowski geometry. 

The numerical accuracy of our scheme can be assessed by comparing the normalization 
of the flux of our calculations for large values of h/M with the analytic result of the flat 
geometry, which for r <C h is F^ at ~ N tot / (Arch 2 ), where N tot is the total number of photons 
emitted by an isotropic source. Our computation proceeds by launching photons starting 
at t/> = and incrementing the launch angle by Aip = 5 x 10~ 5 radians. The total number 
of orbits we compute is then A" = ir/Aip; since each photon is weighted by the sine of the 
emission angle simp (to take into account the number of photons emitted in all azimuths) 
the corresponding total number of photons emitted by an isotropic source, N tot , is related 
to the number of orbits computed, N, by the relation 

N N / \ 2 

N tot = sin(^) = £ sin ^ = cot (^) ~ (for N » 1). (8) 

i=o i=o ^ ' y 

Using this value for A^ot and assuming h/M = 100 (to provide adequate approximation to 
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Fig. 3. — Photon number illumination laws as a function of the normalized disk radius r/M 
for a Schwarzschild and a Kerr black hole with the values of a noted in the figure. In both 
cases the source height has been h/M = 3. Radius of the ISCO is denoted by vertical dashed 
line for a/M = and solid line for a/M = 0.99. The illumination in each case has been 
computed to r = r H . 



Minkowski geometry) one obtains for the flux at M < r < li (using Aip = 5 x 10~ 5 ) the 
value -Ffl at = l/n, in excellent agreement with the value shown in Figure H] for the specific 
value of h. Our scheme therefore, provides the correct shape and absolute normalization 
of the illumination at the appropriate limits and we are confident that it provides also the 
absolute normalization and illumination shape in the regime in which the focusing effects of 
the geometry are important. 

The numerically obtained illumination profiles of Figure H] were fit with an empirical 
function of the following form (in the parameter interval noted) 

Fir h-K) = Kxf Fl(r ' h) f ° r 6 < h/M ~ 100 ' r/M > 115 (9) 
V ' ' ; \ F 2 {r,h) for 3<V M <6, r/M > 1.15 V; 

where K is the normalization and the functions Fi ^{r, h) are given by 

1.026 x 10- 5 (h - 152.9)(/i 2 - 158.7/i + 6569) 



Ft(r,h) 



r 



3 



3.364 x 10- 2 (h - 348.7)(/i + 11.98)(/i + 139.4) 



{h 2 + (l + l/hfr 2 } 



3/2 



(10) 



79.56(/i- 6.250)(/i 2 - 9.763/i + 25.19) 
F 2 (r,h) = 
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Fig. 4. — Radial illumination profiles of the photon number flux at the disk plane for an 
isotropic flux of 40,000 photons at the source. We take a/M = 0.99. The corresponding 
source heights are h/M = 100,50,30,10,6,5,4 and 3, from bottom to top at the left side 
of the figure. Dashed line denotes a slope of —3 for power-law form (r~ 3 ). Plunging region 
inside ISCO is shown by a shaded box. 

3921(/i - 7.364)(/i 2 - 8.556ft, + 21.92) 
{/i 2 + (l + l//i) 6 r 2 } 3/2 

where h and r are measured in units of M; the numerical constants that provide the normal- 
ization of each of the fractions of the fitting formulae were chosen to give the correct normal- 
ization for a total number of photons (by an isotropic source) equal to N tot = 2/ Aip = 40, 000, 
i.e. the number appropriate to our own simulations. 

In Figure the fits provided by this formula (gray) along with the numerically computed 
profiles (solid dark) and the corresponding Minkowski geometry profile (dotted curves) are 
shown for (a) h/M = 100, (b) h/M = 10, (c) h/M = 6, and (d) h/M = 3. It is apparent 
that the analytic formulae provide very good fits to the computed curves both in shape and 
normalization, especially at r < h, the region crucial for the correct determination of the 
relativistic iron line shapes. 

It should be noted at this point the above calculations provide the fraction of the total 
number of photons (number flux) emitted by the source that are intercepted between r and 
r + dr on the disk as measured in Boyer-Lindquist coordinates (t, r, 8, <fi) by an observer at 
rest with respect to the photon source. However, conventional calculations of Fe line profile 
involve the rate of Fe line photon production as a function of the coordinate r, measured in 
the same coordinates by an observer in relative motion with respect to the photon source 
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Fig. 5. — The same illumination profiles (solid dark curves) as in Figure H]for (a) h/M=100, 
(b) 10, (c) 6, and (d) 3 along with our best-fitted analytic function (solid gray curves) 
expressed by equation (J9j). Dotted curves represent the corresponding profiles in Minkowski 
space. 



with the local Keplerian velocity. Modulo an efficiency factor for the conversion of continuum 
into line photons, the transformation of the rates between these two Lorentz frames involves 
powers of the factor 

V=^, (12) 

UeP^ 

where u%, are respectively the four-velocities of the emitter and receiver of radiation and 
p M the photon four-momenta. The expression of this factor for a source at rest in the Boyer- 
Lindquist frame and a receiver at rest with the frame moving with the local (disk) Keplerian 
velocity is computed in Appendix B. The analytic expressions given by equation (jUJ) can 
be transformed to this frame by multiplying by T> 1+a where a is the spectral index of the 
power-law X-ray source (assumed to have the form F v cx v~ a in the local source frame): 
one power of T> indicates the change in the rate photons are intercepted by the comoving 
observer in the fluid/disk frame, while the factor T> a accounts for the spectral shift in the 
photon energy as perceived by the same comoving observer. 
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Fig. 6. — Radial illumination profiles, corresponding to Figure HI of the photon energy flux 
in the comoving fluid frame. Various slopes are shown as a reference. 

In Figure [6] we present the results of this section modified by the above factor (photon 
energy flux) for a Schwarzschild and a Kerr (a/M = 0.99) geometry; i.e. F(r, h; K)T> 1+a . 
We assume a Keplerian motion outside the ISCO while a free-fall motion in the plunging 
region. One sees that in both cases there is a considerable increase in the line production rate 
for sufficiently small radii. Since the slope of illumination in photon energy flux is always 
steeper than that of the photon number flux by the factor of V 1+a , the illumination profile 
of the photon energy flux exhibit two characteristics: the slope of the outer disk zone follows 
that of flat space (q ~ 3), while in the inner disk zone (outside the ISCO) the slope becomes 
(somewhat) steeper (3 < q < 4). Inside the ISCO the slope becomes extremely steep simply 
because of the factor T> 1+a . Assuming that the disks terminate at the ISCO and that there 
is little emission interior to this radius (as shown in Fig. 3 of Reynolds & Begelman 1997 the 
matter is expected to be highly ionized in this region, except in cases of very low efficiency), 
our results indicate that in the case of a Schwarzschild geometry the amplification of the 
photon flux due to the position and motion of the absorbing plasma on the disk makes only 
a small change to the results of Figure HI Howerver, the situation is different in the Kerr 
geometry examined(a/M = 0.99), where the above effect becomes significant at regions of 
the disk outside the ISCO . 

A different aspect of the illumination problem concerns the fraction of the photons of 
an isotropic source that is intercepted by the accretion disk plane. For a point-like source 
at height h on the symmetry axis in Minkowski geometry, the total number of photons 
intercepted in the region rn < r < R of the disk is given by 



N tot - \cos~ 1 (h/VWTB?) - cos-\h/VWTl) 



(13) 



This is the fraction of the total number of photons, iV tot , emitted within the solid angle 
subtended by the disk (excluding the black hole horizon, taken in this case to have the value 
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Fig. 7. — The fractions of the photons emitted by an isotropic source that reach the disk 
plane /p, escape to infinity f^, and lost through the horizon / B h, as a function of h. We 
consider a/M = (solid) and a/M = 0.99 (dotted curves) cases. The cut-off radius is chosen 
to be R/M = 2000. 




t-r/M = 1.15 ~ 1). For R ^> h ^> 1 this attains the value N tot /2, the remaining photons 
escaping to infinity. The focusing effects of the geometry modify significantly the fraction 
of photons intercepted by the disk, escape to infinity or are absorbed by the black hole 
as the source approaches the horizon. The procedure we have used in §2 can be used to 
compute these fractions which are of value in estimating the effects of the geometry on the 
ratio of line flux to continuum, i.e., the line equivalent width (EW) or the strength of the 
Compton reflection feature. Figure [7] shows the percentage of the photons of an isotropic 
source that reach the disk plane (out to a distance R/M = 2000), /d, that of photons 
escaping to infinity, f^, as well as the fraction absorbed by the black hole, /eh, as a function 
of the source height h for a/M = (solid) and a/M = 0.99 (dotted curves). For large h 
{h/M > 50) /d, foo — > 1/2 while /bh — ► 0, as expected. As the height of the source decreases, 
foo decreases monotonically to zero (for h — > th); however, /□ increases initially reaching a 
maximum of ~ 70% for h/M ~ 3 (~ 60% at h/M = 5 for a Schwarzschild black hole) before 
it declines to zero as the source approaches the black hole horizon. Last, the fraction /bh 
does have the qualitative behavior expected, being negligible for large h, reaching unity as 
h — > rn- The qualitative behavior of these quantities does not change as one moves from 
the case of an extreme Kerr geometry (dotted) to that of a Schwarzschild black hole (solid 
curves). However, the maximum value of /p is slightly smaller and it is achieved for a slightly 
larger value of h. It i s of interest to note that the values of /bh, /d, foo are very similar 
to those presented by iMiniutti et al.l (120031 ). who considered a different source geometry (a 
ring-like source in an extreme Kerr geometry of radius r/M = 2 located at a height h) 
and a different computational approach. It appears therefore, that this represents a generic 
behavior for illumination not only by point-like but also by extended X-ray sources provided 
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that they are concentrated near the symmetry axis of the configuration (rather than, e.g., 
near the equator). 



3.2. Illumination by an Anisotropic X-ray Source 

While the formulae given above provide very good approximations to disk illumination 
by an isotropic source, we have chosen to produce also the illumination by an anisotropic 
emission of photons, appropriate for cases that the observations demand illumination more 
concentrated near the disk inner edge than that of an isotropic source. Clearly this is an 
arbitrary choice for if one moves away from the isotropic assumption there is little constraint 
as to how much anisotropy should be used. Nonetheless, in order to provide an algorithm 
as well as a gauge of the source anisotropy with a certain theoretical basis we have decided 
to Lorentz-boost a locally isotropic source distribution along the symmetry axis and toward 
the black hole by a constant velocity v. We assume the emission to be instantaneous, i.e., 
that the photons are emitted before the source has moved substantially to alter its height h. 

The methodology followed is similar to that used earlier: a value of the local angle ip' 
is chosen uniformly as before to provide the photon distribution at the rest frame of the 
emitter. We then compute the corresponding angle if) in the local rest frame which is related 
to ip' by the usual angle transformation 

sin?// 

tam/> = -T— - , 14 

7(cos-7// — p) 



where (3 = v/c and 7 = 1/y 1 — (3 2 . This value of the angle is then used along with equation 
([6]) to obtain the values of b and f used in the integration of the photon orbits. The factor 
sin^ of this latter angle ip is again used to provide the correct normalization of the photons 
illuminating the disk. 

As expected, the anisotropy of emission, measured by the velocity of the corresponding 
relativistic motion, provides a higher concentration of the available photons towards the 
inner edge of the disk than the geometry alone for a source approaching the disk ((3 > 0) 
and correspondingly less concentrated for a source moving away from the disk (/3 < 0). At 
the same time, depending on the anisotropy and the height of the source, a large fraction of 
the photons can be lost through the black hole horizon. 

Figure [8] shows the photon trajectories by an anisotropic source (moving toward the 
disk) of anisotropy corresponding to instantaneous velocities (a) (3 = 0.5 and (b) (3 = 0.95 at 
a height h/M = 6 illuminating a disk around a black hole of a/M = 0.99. Figure [9] displays 
the corresponding disk illumination of an anisotropic source for two different values of the 
source height, namely (a) h/M = 6 and (b) h/M = 10. For each value of h/M we consider 
for comparison source anisotropies corresponding to (3 — —0.95, —0.5, 0, 0.5, 0.95, with 
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(a) p = 0.5 
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Fig. 8. — The photon trajectories in the poloidal plane for an anisotropic source at h/M — 6 
with (a) (3 = 0.5 and (b) 0.95. 



the minus sign indicating motion away from the accretion disk. As expected, the anisotropy 
increases disk illumination for r < h for sources approaching the disk, while decreasing it 
for sources receding from the disk; its slope becomes correspondingly steeper for the higher 
value of the anisotropy parameter (3 = v/c. As a result, the profile can deviate greatly from 
that of an isotropic source. Qualitatively, the increase in (3 mimics a decrease in h to the 
point that cos ift ~ I/7 = a/1 — (3 2 , that is to the point that the Lorentz boosting solid angle 
(oc I/7 2 ) is roughly equal to the solid angle subtended at the source position by the black 
hole horizon; at this value of f3 the illumination can be approximated by a power-law r -3 for 
the entire r— range as shown in Figure [9J correspondingly a negative value of f3 (a receding 
source) provides an illumination profile equivalent to a source at a higher value of h. We have 
checked by direct computation that further increase in the anisotropy through an increase in 
(3 (or 7) preserves the r -3 functional form of the illumination (at all distances) but decreases 
its normalization in proportion to 7 -2 as might be expected on the basis of qualitative 
considerations; by contrast increasing of \(3\ for a receding source, increases indefinitely the 
effective height of the source h. Cases that demand an illumination law steeper than r -3 (in 
number flux) must therefore involve assumptions about the source and its geometry different 
from those hitherto considered (see §4). On the other hand, energy flux profile can become 
steeper than ~ r~ 3 by the factor of T> 1+a in energy flux, as explained earlier. In this respect 
it is worth considering that a source that could reverse its axial motion (±/3) would lead to 
substantial change in the number of photons intercepting the accretion disk and the flux of 
the Fe line, even for moderate values of \/3\ (~ 0.5). 
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Fig. 9. — The illumination profiles for (3 = —0.95 (gray), —0.5 (gray), 0,0.5,0.95 from 
bottom to top at the left side of the figure for (a) h/M = 6 and (b) 10. Three lines represent 
the slopes for —2 (dotted), —3 (solid) and —4 (dashed lines) as a reference. 




Fig. 10. — Same as Fig. [7]but for anisotropic cases. A source approaching the black hole has 
P > 0. The source height is set to be h/M = 6. 



For the reasons discussed in the previous subsection we have also computed for the 
anisotropic illumination case the fractions /bh, Id, foo of the photons intercepted by the 
black hole, the accretion disk and escaping to infinity, applying the same methodology used in 
§3.1. These results are shown in Figure [TOl for a source at a fixed height h/M = 6 as a function 
of the anisotropy parameter (3 = v/c for a/M = (solid) and a/M = 0.99 (dotted curves). 
For p = these fractions are those corresponding to h/M = 6 of Figure [TJ as they should. As 
shown in this figure the maximum (integrated) disk illumination, /p, is obtained for P ~ 0.3. 
Also as expected, /» decreases with increasing P while /bh increases approaching unity as 
P — > 1. These fractions are very different from those of flat space, which can be calculated 
analytically by integrating the appropriate photon angular distribution [oc 7~ 2 (1 — Pn)~ 2 } 
over the cosine \x of the polar angle 6 from fi = — 1 to \i = (for foo), from \i — to 
/i c = h/y/1 + h 2 (for /d) and from /x c to /x = 1 (for /bh) where /i c is the cosine of the angle 
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that the horizon subtends at the distance of the source. In the flat geometry one expects the 
maximum disk illumination to be obtained for 9 ~ 3 cos~ 1 (/x c ) or (3 ~ cos[3 cos _1 (// c )] ~ 0.85, 
a value substantially different from that given in Figure (TUl the difference attributed to the 
focusing effects of the geometry. Also, the ratio of /d to f^, a quantity relevant for the EW 
of a line, achieves in flat space its maximum value for (3 ~ 0.65. 



3.3. The Iron Line Profiles 

In this section we produce the relativistic Fe line profiles corresponding to the illu- 
mination laws derived above. The pr ocedure followed for that is the standard ray-tracing 
method given in the li t erature (e.g., Fanton et al. 1997 ; Cadez. Fanton. fe Calvanil Il998l ; 



Beckwith fe Pond 12004 ; iFukumura fe Tsurutal 12004 ). by which the photon orbits are fol- 



lowed from the source to the observer's field of view (or image plane). In all cases the 
inclination angle of the observer is 30°, relevant for Seyfert 1 galaxies, and the energy resolu- 
tion is assumed to be AE = 150 eV. We only consider regions outside the ISCO for producing 
line photons. In Figure fTTl sample line spectra are shown for (a) isotropic and (b) anisotropic 
sources. In (a) the line profiles are shown for h/M = 3, 6 and 10 to illustrate the effects of 
changing the source height on the line profile for (3 = 0. In (b) we show the line profiles for a 
source at a fixed height h/M = 6 but of varying anisotropy given by (3 = —0.95, —0.5, 0, 0.5 
and 0.95 to exhibit the effects of source anisotropy on the corresponding line profiles. 

In (a) the profiles are generally broader the closer the source is to the disk due to the 
increased focusing of the photons toward the disk inner edge. Because the line profiles were 
produced by collecting photons only out to disk radius r out /M = 30, the line flux decreases 
with decreasing h/M; this is because for a decrease in h/M, the photon flux at r/M > 10, 
which contributes to the narrow peak of the line near 6.6 keV, decreases, while it increases 
for the part of the disk that contributes to the low (E < 5 keV) emission. At the same time, 
however, as shown in Figure [7] the fraction of the photons escaping to infinity also decreases 
and as a result the line EW will not necessarily decrease (provided that the point source is 
the only source of illumination). 

In (b) the calculated line spectra become broader with increasing anisotropy of an 
approaching source (0 < (3 < 1), simply because more photons are concentrated near the 
inner edge of the disk. For sufficiently high values of (3 > 0, the anisotropy beams a large 
fraction of the photons towards the black hole horizon, resulting in a decrease of the line 
absolute flux. The red wing of the line, on the other hand, is not suppressed with non-zero 
flux at energies less than 4 keV. However, as discussed above, the number of photons escaping 
to infinity (the observer) also decreases and so the EW of the line may not necessarily be 
negligible. Figure [Tib exhibits also, as a reference, the line profile in the same geometry 
produced by sources receding from the disk (/3 < 0). The beaming effect in this case occurs 
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Fig. 11. — Iron line spectra (in arbitrary units) for (a) isotropic and (b) anisotropic sources 
with a/M = 0.99. In (a) we show h/M = 3, 6 and 10 for (3 = while in (b) we use 
(3 = —0.95, —0.5,0,0.5 and 0.95 for h/M = 6. Vertical dashed line denotes a cold iron line 
at 6.4 keV in the rest frame. 

in the opposite direction (away from the black hole) suppressing the overall line flux at all 
energies. 

To provide an estimate of the influence of the source position and anisotropy on the line 
EW we present in Figure [T21 the ratio of the photons escaping to infinity to those intercepting 
the disk /oo//d as a function of (a) the source height h/M for the case of an isotropic source 
and (b) the anisotropy parameter j3 for the case of an anisotropic source at h/M = 6, i.e. 
cases that correspond to Figures [7] and [101 respectively. While both fluxes are reduced 
as the source approaches the black hole (or its anisotropy increases) the ratio of the fluxes 
intercepted by the disk and escaping to infinity actually increases. All other things remaining 
equal (e.g. efficiency of Fe line production), this should lead to an increase in the line EW. 



4. Discussion & Conclusions 



We have produced through direct integration of a large number of photon orbits the 
illumination law of thin accretion disks around black holes by point-like sources at different 
heights above the black hole on the axis of symmetry. We have also produced analytic fitting 
formulae that can provide accurate approximations to the illumination laws (in number flux) 
under the above assumptions for both Schwarzschild and Kerr black holes. The corresponding 
energy flux illumination profiles can be obtained by multiplying those by the enhancement 
factor, T> 1+a approriate for each situation considered. These can be employed as needed to 
calculate the corresponding line profiles in a fashion similar to that shown in §3.3. We have 
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h/M p 



Fig. 12. — The ratio of the fractions of the source photons escaping to infinity to that 
intercepted by the disk /d for (a) an isotropic source (given as a function of the source 
height h/M) and (b) for an anisotropic source (given as a function of the anisotropy velocity 
(3). Panel (a) corresponds to the flux ratio of Figure [7] while (b) corresponds to the flux ratio 
of Figure [TU1 Solid lines are for a/M = while dotted ones for a/M = 0.99. 

also produced for each case the corresponding fractions of photons emitted by the source 
that escape to infinity, intercept the disk or being absorbed by the black hole, quantities 
necessary in producing the relative normalization between the escaping (observed) flux and 
that reprocessed into either the Fe line or the Compton reflection feature. Our overall results 
can be summarized as follows: 

(a) For sufficiently large values of h/M (> 50) the illumination profiles of isotropic 
sources approach those of flat space, with the focusing effects of the geometry becoming 
noticeable for r/M < 7. As the height h approaches the size of the event horizon, the 
focusing effects of the geometry concentrates the photon paths to the inner edge of the 
accretion disk while reducing the illumination of larger radii; these features bear direct 
relation to the shapes of the resulting line profiles. 

(b) In addition to the focusing effects of the geometry, of interest in the production of Fe 
line photons by the incident radiation is also the increased rate of photon reception and the 
photon energy shift between the emitting source and the receiving plasma, on a Keplerian 
orbit on the accretion disk. Assuming that the disk and line production is important only for 
radii larger than that of ISCO, this amplification is significant for an extreme (a/M = 0.99) 
Kerr geometry but not for a Schwarzschild one. As shown in Figures H] and [6] the main 
effect of the hole spin parameter a is to reduce the size of the horizon and the size of ISCO, 
which simply determines the minimum radius at which the illumination (and hence the 
contribution to the line profile) needs to be calculated. This is a rather surprising outcome 
of our calculations; despite the different topologies of the orbits with increasing black hole 
spin a the resulting illumination laws appear to change very little for h/M > 3 (however the 
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profiles become significantly different for smaller values of h). As a result the fitting formulae 
we provide in §3.1 (with the modifications for the rate and photon energy discussed in §3.1) 
should suffice to produce illumination profiles by isotropic sources for both Schwarzschild and 
Kerr black holes. 

(c) Motivated by the need of illumination more concentrated to the disk inner radius, in 
order to explain certain particularly broad Fe lines, we have modeled the illumination of disks 
by point-like anisotropic sources. We modeled the source anisotropy by Lorentz-boosting an 
isotropic source along the symmetry axis and toward the black hole by a constant velocity 
p. We found that for sufficiently large values of (3 > 0, the illumination law (in number flux) 
assumes the shape of a single power law oc r~ 3 , over the entire range of radii. We have also 
found that further increase in /3 does not modify this functional form but simply decreases 
the overall normalization because of absorption of most of the photons by the black hole. 
While the black hole spin does not affect the illumination profile it does affect the ratio of 
source photons intercepted by the disk (which can contribute to the line emission) to those 
that escape to infinity and contribute to the continuum emission, in particular in the case 
of anisotropic source illumination. 

(d) We found that the photon energy flux illumination law, obtained by multiplication 
of the photon number flux by the factor V 1+a , can be significantly steeper in the innermost 
disk region of a Kerr geometry with an index that approach the value q ~ 4 for a small range 
in radius. 

Our calculations indicate that it is hard to obtain disk illumination profiles more concen- 
trated than those discussed above (demanded by certain observations) within the confines 
of our assumptions. However, modification of these assumptions should easily allow for 
much steeper profiles. A point-like source on the symmetry axis is an idealization that 
simplifies the computational task. A point-like source off axis will greatly enhance the illu- 
mination in the region immediately underneath its position; this requires a more complicated 
non-axisymmetric description, not demanded as yet by the observations (but it may in the 
future). It is also apparent that an extended source can produce arbitrarily steep illumina- 
tion profiles depending on the assumption of the radial distribution of the X-ray intensity 
within the source. In the absence of any physically compelling prescription, it is simplest to 
approximate the illumination as a combination of power-laws: (a) outer illumination oc r~ 3 
appropriate for large distances; our results indicate this to be a valid assumption for dis- 
tances much larger than the vertical extent of the source and (b) a steeper one appropriate 
for the illumination of the inner portions of the disk. While the computation of the resulting 
illumination is possible under an extended source assumption by extending the procedure 
described in this note we do not believe that such a model would lead to any significant 
insights about the physics of these sources. 

To this point we have not discussed at all the nature of the illuminating source. Several 
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proposals have been made in the literature of w hich we provide an incomplete list: An accre- 
tion disk corona (e.g.. lHaardt &: Maraschilll99ll ). powered either by magnetic fields threadin g 



the disk or by magnetic fields threading the black hole (e.g., lLill2002l ; iReynolds et al.l 12004 ) . 
An alternative is the emission by shocks in the accre tion flows that funnels (through magnetic 
fields) the disk material onto the black hole (e.g., iFukumura. Takahashi. fe Tsurutal 120071 . 
for magnetohydrodynamic shocks in accretion). The source can be either sta tionary or vari 
able, as such shocks w ere proposed to be unstable and oscillating (see, e.g.. iTrussoni et al. 
19881 ; lAoki et al.l 12004 for discussion of shock stability). Finally, we have not addressed 
at all the issue of disk ionization, an issue that affects its line produc ing efficiency (e.g., 
Nayakshin. Kazanas. fe Kallmanll2000l : iBallantyne. Ross, fe Fabianll200ll ) as they are not re- 
lated directly to the main thrust of the present work. 

Assuming that the plunging region (inside the ISCO) is optically-thin due to very low 
mass-accretion rate, photons emitted from the source could in principle cross the plunging 
region. We have followed photon trajectories (from the source to the equatorial plane) and 
found that certain of them turned around (due to strong light bending) to hit the bottom side 
of the disk at larger radii and contribute to the illumination there; we found this contribution 
to the illumination to be unimportant. 

We have produced in §3 fitting formulae for the illumination of accretion disks (in 
photon number flux), which can be easily converted into the photon energy flux needed in 
the calculations of Fe line profiles. Our models and fits are consistent with the values implied 
by spectral fits to the Fe lines of most AGNs. Assuming power-law distributions oc r~ q for 
emissivity, the values of the index q obtained fo r MCG-5-23-13, NGC 4051 and NGC 3516 in 
recent Suzaku observations (IReeves et al.l 120061 ) are q ~ 2 — 3; our models can accommodate 
easily the lower values of q while demanding a reasonably high degree of anisotropy for values 
q > 3. However, the disk illumination under the assumptions considered in this note cannot 
account for the width of the Fe line profile of MCG -6-30-15 as obtained ei ther by XMM- 



Newt on (jWilms et al.ll200ll ; iFabian et al.l 120021 1 or more recently by Suzaku (jMiniutti et al. 
20071 1 which demand a very steep (q ~ 4.6) illumination of the inner disk section along with 
a flatter one (q ~ 2.6) at larger radii. As we discussed earlier, the easiest scheme that could 
accommodate such a profile would be an extended (r/M ~ 2) rather than the point-like 



source used in our models. An example of such a source has been given in iMiniutti et al. 



( 120031 ). however we believe that other types of extended sources could provide good fits as 
well. Another possibility has also been proposed that m agnet ically-induced torque (on the 
disk) could produce very steep emissivity profiles (e.g.,|Lj|2002|). Modeling such sources goes 
beyond the scope of the present work and maybe undertaken in the future. 
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5. Appendix A 

In the Appendix A we present the orbit equations which we have integrated in producing 
Figures [2] and [S] in §3. The symmetry of the problem defines an axis which in general is 
different from the axis connecting the source to the black hole. In this case, in addition to 
the total angular momentum £ of the photon with respect to the black hole, there is also 
its component about the hole's rotation axis L z . Therefore, in this case the photon has two 
degrees of freedom, i.e. the impact parameters with respect to the black hole rj = C/E 2 and 
the impact parameter with respect to the hole's rotation axis £ = L z /E. 

In this case the equations of motion in the r and 9 directions are 

and * = ^>, (15) 



This preprint was prepared with the AAS IATgX macros v5.2. 
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where E = r 2 + a 2 cos 2 8 and A = r 2 — 2Mr + a 2 , with a being the angular momentum per 
unit mass associated with the Kerr metric, and the functions R(r) and 0(9) are defined by 

R(r) = \(r 2 + a 2 )E - aL z } 2 - A[C + (L z - aE) 2 } , (16) 
0(9) = C - (-a 2 E 2 + L 2 cosec 2 #) cos 2 fl . (17) 



In analogy with equation (j2j) these are supplemented by the following equations for t and 

■ _ dt _ (r 2 + a 2 ) [(r 2 + a 2 )E - aL z ] + aAL z - a 2 EA sin 2 # 
' = ~cfr~ A(2r 2 + a 2 + a 2 cos29) ' 



• _ d(p _ a[a 2 E - aL z + E(r 2 - A)] + L Z A cosec 2 6> 
^ = rf7 = 2 A(2r 2 + a 2 + a 2 cos29) ' ^ 

By dividing equations (fT6"j) and (|17[) above by E 2 and re-expressing them in terms of the 
impact parameters r\ and £ we obtain 

R(r) = l(r 2 + a 2 )-at} 2 -Al V + (Z-a) 2 }, (20) 
0(9) = 7] - (-a 2 + £ 2 cosec 2 #) cos 2 # . (21) 

Specializing to a source located on the axis of symmetry, £ = 0; then the equations of 
motion and the integrals of motion read respectively as 

S 2 r 2 = (r 2 + a 2 ) 2 - A(rj + a 2 ) , (22) 
Y?9 2 = v + a 2 cos 2 9 , (23) 

and 

(r 2 + a 2 ) 2 -a 2 Asin 2 fl 

A(2r 2 + a 2 + a 2 cos2#) ' 1 ' 

■ a(a 2 + r 2 - A) 

9 A(2r 2 + a 2 + a 2 cos2#) ' 1 ' 

The angle between the axis perpendicular to the disk and the direction of photon emission 
(we call it 5 instead of ip in this case) is given by an expression equivalent to that of equation 

\l/2+ 1 



cot 5 



\9 



rr 



11/20 A 1 / 2 



(r 2 + a 2 ) 2 -A(^ + a 2 )] 1/2 

(26) 



1] + a 2 cos 2 9 

which can be readily inverted to provide an expression for the impact parameter r\ in terms 
of S, namely 

v (5, 9; r,a) = -^ [(r 2 + a 2 ) 2 - a 2 A(l + cos 2 # cot 2 5)] sin 2 5 , (27) 
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an equation that reduces to equation (jfJJ) for a = 0. 



We follow a procedure similar to that of the Schwarzschild case, i.e. we opt for the 
numerical integration of second-order equations for r and 9. Being emitted from a point 
along the symmetry axis (£ = 0) the photons still have one degree of freedom namely the 
impact parameter 77. We use equation (1271) to relate the local emission angle 5 to rj and 
equations (j2"2"l) and (1231) to obtain the initial values of r and 9 necessary for the integrations. 
The integrations proceed again by a choice of the angle 5 chosen uniformly between and 7r 
from (r = h, 9 = 0) until 9 = tt/2 and r is greater than the black hole horizon corresponding 
to the specific value of the hole angular momentum a. For completeness we present below 
the second-order equations used in the integrations. These are the following: 



1 

S2 



(a 2 sin 29 9 - 2 rr)rS + 2(r 2 + a 2 ) r - (r - M) (r) + a 2 ) 



and 



9 



1 



2S 2 



2(a 2 sin 299 - 2 rr)9E - a 2 sin 29 



(2? 



(29) 



The general case for a photon source off the axis, £ 7^ 0, is straightforward to obtain using 
the general form of the equations for R(r) and 0(9). However, in this case one needs to 
choose an additional angle in the (local to the source) azimuthal direction which has to be 
related to £ in a similar fashion. In that case the illumination pattern depends also on the 
azimuthal angle <fr and the results sufficiently complicated to be presented within this work. 



6. Appendix B 

As emitted photons hit the accretion disk, its local energy and the rate at which they 
are intercepted by the disk will be relativistically shifted. In order to calculate photon energy 
flux measured in fluid/disk comoving frame, one must take those effects into account. Photon 
redshift factor from an emitter's frame to a receiver's frame is defined as 



V 



v u 'P) 



UeP^ 



(30) 



where "r" 



and "e" respectively denote the receiver and emitter. For photon's four-momentum, 
we hav e p = (— E, ±E\/R(r )/ A, ±y / 0(^), £E) where E is the photon energy in the emitters 
frame (jChandrasekharlll983l . p. 347). For four- velocity we generally assume u = (u*, u r , 0, u^). 
Thus, we obtain 



X> = I -7 



u' 



;i-(j] T yfj//A) r 

:i-& T vw)v r /A) e 



(31) 
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where v r = u r ju l and Q = u^/u 1 . For example, for photons emitted from a stationary 
on-axis source (£ = 0) reaching the accretion disk (outside ISCO), we find 



V 



u 



disk 



(g tt + 2g t <fl K + g<f>4>Q x) disk 



1/2 



(32) 



where Qk 
p.336). 



M/(avM + r 3 / 2 ) is Keplerian frequency of the disk ( Chandrasekharl Il983 



